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■The Texas A&M University Penning Trap facihty is an upcoming ion trap that wiU be used to search for possible scalar 
^currents in T = 2 superallowed /3-delayed proton decays, which, if found, would be an indication of physics beyond the 
■standard model. In addition, TAMUTRAP will provide a low-energy, point-like source of ions for various other applications 
^at the Cyclotron Institute. The experiment is centered around a unique, compensated cylindrical Penning trap that 
■employs a specially optimized length/radius ratio in the electrode structure that is not used by any other facility. This 
[allows the geometry to exhibit an unprecedented 90 mm free radius, which is larger than in any existing trap, while 
■at the same time remaining a tractable overall length. The trap geometry was designed from first principles to be 
[suitable for a wide range of nuclear physics experiments. In particular, the electrode structure is both "tunable" and 
■ "orthogonalized" , which allows for a near quadrupole electric field at the trap center, a feature necessary for performing 
[precision mass measurements. 
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1. Introduction 

Low energy precision /3-decay experiments have proven 
.to be an excellent compliment to high energy physics for 
'placing new constraints on physics beyond the standard 
[model (SM) Up to this point, it has been possible 

"to explain the results from such experiments by a time 
.reversal invariant V — A interaction displaying a maximal 
■violation of parity; however, more precise measurements of 
.ft values [1,1^ and correlation parameters Q in particular 
■/3-decays can serve to test for the presence and properties 
.of any non-SM processes that may occurr in such interac- 
■tions. 

1.1. Motivation 

The initial experimental program at the upcoming Texas 
A&M University Penning Trap (TAMUTRAP) facility will seek 
to improve the limits on non-SM processes in the weak in- 
teraction, in particular scalar currents, by measuring the 
P — V correlation parameter, ap^^ for T = 2, 0+ ^ 0+ 
superallowed /3-delayed proton emitters (the generic decay 
scheme is shown in Fig.[Tl and the preliminary list of nuclei 
to be studied is outlined in Table [T]). 
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Figure 1: Decay scheme for a generic T : 
/3-delayed proton decay. 




: 2, 0+ — > 0+ superallowed 



The general /3-decay rate with no net polarization or 
alignment is given by [Tf: 



dEpdflpdQ, 



oc 1 + — a^i, cos( 



(1) 



where, E^, Pe, and nie are the energy, momentum, and 
mass of the /3, Op^, is the angle between the /3 and j/, and 
b is the Fierz interference coefficient. Thus, it is possible 
to determine the (3 — v correlation parameter by means of 
an experimental measurement of the angular distribution 
between the (5 and v. For the strict V — A interaction 
currently predicted by the SM (in pure Fermi decays) the 
j3 — V angular distribution should yield a value for a^n 
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24Si 


147.15 


3.91 


40.8 


28 g 


180.33 


3.70 


39.7 


32 Ar 


141.38 


3.36 


37.8 


36Ca 


141.15 


2.55 


33.0 


40Ti 


72.13 


3.73 


39.9 


48Fe 


63.48 


1.23 


22.9 



Table 1: The T = 2 nuclei that will compose the initial experimental 
program measuring a^,^. The Larmour radii, Ri^, for the ejected 
protons of interest (having energy Ep) shown are calculated for the 
7 T magnetic field of TAMUTRAP. 



of exactly 1. Any admixture of a scalar current to the 
predicted interaction, a result of particles other than the 
expected being exchanged during the decay, would 

result < 1. 

TAMUTRAP will observe this angular distribution between 
f3 and v for /3-delayed proton emitters in order to take 
advantage of the benefit these particular decays have on 
the experimental procedure. In such a case, the /3-decay 
yields a daughter nucleus that is unstable, and can result 
in the subsequent emission of a proton with significant 
probability. As discussed in 6], the great advantage to 
utilizing /3-delayed proton emitters for such a study is that 
this proton energy distribution contains information about 
0pi,. If the /? and v are ejected from the parent nucleus in 
the same direction, they will impart a larger momentum 
kick to the daughter nucleus, which will be inherited by the 
proton. Conversely, if the /3 and v are emitted in opposite 
directions, this momentum kick is reduced. By measuring 
the proton energy distribution at TAMUTRAP the value of a/ji^ 
will be deduced, which can then indicate the existence of 
scalar currents in these decays 01. 

1.2. Cylindrical Traps 

A cylindrical Penning trap geometry allows for efficient 
access to the trapped ions, a large trapping volume com- 
pared to a hyperboloid trap geometry, and an electric field 
that can be described analytically, which is of particu- 
lar importance during the design process. Additionally, 
cylindrical electrodes are more easily manufactured with 
higher precision. For these reasons, Penning traps with a 
cylindrical geometry have been widely employed in nuclear 
physics research experiments ranging from precision mass 
measurements [1| to the production of anti- hydrogen . 

Precision /3-decay experiments are well served by a 
cylindrical geometry due to the fact that the magnetic field 
employed to trap the ions radially may simultaneously be 
used to contain charged decay products [13] such as /3's 
and protons with up to 47r acceptance in an appropriately 
designed trap. The strong magnetic radial confinement in 
combination with the weak electrostatic axial confinement 
direct the decay products of interest to either end of the 
trap for detection with negligible affect on the energy of 



the particles. At the same time, features of a cylindrical 
trap geometry can be useful for other nuclear physics ex- 
periments, such as maintaining a line of sight to the trap 
center for spectroscopy, an easily "tunable" and "orthogo- 
nalized" electric field for experiments requiring a harmonic 
potential (such as mass measurements), and unrivaled ac- 
cess to the trapped ions. 

2. Design of TAMUTRAP 

For the reasons mentioned in ijl.21 a cylindrical geome- 
try has been chosen for the TAMUTRAP measurement Penning 
trap Uj . This particular geometry has been optimized to 
create a design that is suitable both for the precision (3- 
decay experiments of interest, as well as a wide range of 
nuclear physics experiments as discussed. Specifically, the 
trap must display a large-bore for containment of decay 
products, it should allow for the placement of biased de- 
tectors at both ends for observation of these products, and 
it should exhibit a "tunable" and "orthogonalized" geom- 
etry in order to achieve a harmonic electric field. 



2.1. Large-bore 

For (S — v correlation measurements, the trap must have 
a free diameter large enough to contain the decay products 
of interest within the electrodes via the Lorentz force im- 
posed by the trapping magnetic field. The initial program 
of measuring a^^ will investigate the T = 2 nuclei shown in 
Table [1] by observing the proton energy distribution. To 
fully contain protons of interest with full An acceptance 
the trap radius is set to twice the Larmour radius of the 
most energetic expected proton within the 7 T magnetic 
field provided by the Agilent 7T 210 ASR magnet ^\. 



The trap radius was chosen to be 90 mm, which will fully 
contain protons of up to 4.75 MeV. This radius will be the 
largest of any existing Penning trap and will easily contain 
the protons of interest, as well as the less magnetically rigid 
/3's, for the initial a^^ studies. 

2.2. Endcaps 

The other primary requirement for performing the men- 
tioned correlation measurements is that the design must 
ultimately accommodate position sensitive strip detectors 
at either end of the trap. The charged decay products (/3, 
p, daughter) exhibit Larmour precession contained com- 
pletely within the bore of the Penning trap until they are 
detected by such a detector. These detectors have been 
simulated separately as disk-shaped "endcap electrodes." 
Such an approximation satisfies the need to bias the de- 
tectors at some arbitrary potential not necessarily equal to 
that of the hollow cylindrical end electrodes. The impact 
of the endcap electrodes on the solution to the complete 
electric field is discussed in detail in 
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Figure 2: The trap geometry and boundary conditions used in tlie 
analytic solution. 



2.3. Tunable 

A priority for the TAMUTRAP facihty is to attain a very 
good quadrupole field at the center of the measurement 
Penning trap. To achieve this, the trap design must be 
"tunable" , that is, it must make use of compensation elec- 
trodes that serve to adjust the field shape. Other cylindri- 
cal traps 13, 13] already employ tunable geometries; how- 
ever, these configurations were not suitable for TAMUTRAP 
due to the large-bore requirement: enlarging such geome- 
tries (which does not inherently affect the field shape if 
all features are scaled appropriately), results in a trap too 
long to fit within the available 7 T magnet. In addition, 
the analytic solutions of the electric field used to design 
these existing traps [l^ directly employ long-endcap ap- 
proximations, which makes a new design utilizing the ex- 
isting field calculations not applicable to the short-endcap 
/ large-bore requirement necessitated by the envisioned 
TAMUTRAP experimental program. 

For these reasons, a new analytic solution for a short- 
endcap, tunable Penning trap was derived from first prin- 
ciples, in part following the discussion in Ref. [l^. The 
contribution to the potential due to the various electrodes 
(Fig. [2]) can be found by noting that any potential may be 
expanded in terms of the Legendre polynomials, Pk, the 
potential depth, radius from the trap axis, and character- 
istic trap distance [15| , Vq, r, and d, respectively, and the 
expansion coefficients, Ck (here k is even due to symmetry 
across the trapping plane): 
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By superposition, the potential at the trap center may 
also be written as a sum of the potentials of each of the 
contributing electrodes: 



V = Vb0O + ^101 + 



(4) 



Here, (/>o is due to the boundary conditions generated by 
the ring and end electrodes (the primary potential well), 
and (j)i and 4)2 are from the compensation and endcap 
(detector) electrodes respectively (Fig. [5]) , while the Vi 's 
are the potentials at which these electrodes are held. Each 
of the 0's can in turn be expanded in Legendre polynomials 
and substituted back into Eq. 01 which yields the simple 
result: 

V^ V) 

Ck = Dk + Ek-^+Fk-^, (5) 

where Dk, Ek, and Fk are the individual expansion coef- 
ficients due to the ring and end electrodes, compensation 
electrodes, and endcap (detector) electrodes respectively, 
and must be solved for individually. and , which are 
due to hollow cylinder shaped electrodes, can be found by 
expanding 0o a-nd 0i in Bessel functions, Jq. After apply- 
ing the periodic boundary condition in z, (pi^) = 4'{~z), 
it can be found that 
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where An is an additional expansion coefficient. Here fc„ 
is due to the periodic boundary condition, and is given by 
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where ztot = Zr + Zc + Ze + 3zg. Setting Eq. [5] equal to the 
expansion in Legendre polynomials (Eq. [5]) allows one to 
solve for the expansion coefficients of the Legendre poly- 
nomials by equating the two along z. This results in the 
following solutions: 
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The coefficients are subsequently determined by ap- 
plying the appropriate boundary conditions (Fig. [2]) along 
with the orthogonality of cosine, yielding 
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Ztot nPo) 
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The contribution to the potential from the endcap elec- 
trodes must be handled differently. Here <f)2 is defined at 
Ztot for any radius less than po- After simplifying due to 
the cylindrical symmetry of the system, (/)2 can be written: 



(10) 



n=0 



where, kmn is related to Xmm the n zero of the m Bessel 
function as in 
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Bn can now be determined through the application of ap- 
propriate boundary conditions (see Fig. [2]) and the Bessel 
function orthogonality relation, giving 
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Taking into account both endcaps (at ±z) yields the com- 
plete formulation of the potential due to the endcaps at 
any z: 
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This result for the potential was subsequently Taylor ex- 
panded using Mathematica J[Q], yielding coefficients which 
may be referred to here by T^. Setting the result from 
the expansion equal to the potential expanded in Legen- 
dre polynomials along the z-axis and equating terms yields 
the final result for the coefficients Fk'- 



(14) 



All expansion coefficients from Eq. [5] have now been com- 
pletely defined, and the electric field at the trap center can 
be specified to arbitrary precision. By construction, it is 
easy to characterize the components of this field, which 
allows for a straightforward optimization. 

2.4- Orthogonal 

For certain experiments (such as precision mass mea- 
surements), it is crucial to be able to tune out the anhar- 
monic terms (C>4 0) of the electric field during the 
course of a measurement without affecting the harmonic 
(C2) component of the field. In order to achieve this for 
TAMUTRAP, the procedure discussed in Ref. [11] was followed. 

The strength of the dominant anharmonic term for the 
superposition of the electric fields is given by the coefficient 
C4, and the second most dominant contribution comes 
from Cg, which has an affect on the shape of the elec- 
tric field smaller than C4 by a magnitude of (r/d)^ (where 
r is that radius of ion motion and d is the characteristic 



trap dimension) [15[. These coefficients are determined, 
in turn, by the anharmonic contributions from the con- 
stituent electrodes, that is, D4, Dg, E4, Eq, F4, and Fq. 
C4 may always be tuned out by adjusting the potential on 
the compensation electrodes until the field is essentially 
harmonic. For a general geometry, however, this proce- 
dure will affect the value of C2. Since only the potentials 
of the compensation electrodes are adjusted in this pro- 
cess, it is possible to eliminate this affect by requiring that 
E2 = 0, that is that the compensation electrodes have no 
influence on the harmonic term of the superposition, C2 ■ 
The expansion coefficient E2 , which is a function of the 
entire geometry, was minimized using the analytic solution 
derived above. To do this, the constraints imposed by 
other considerations (trap bore, assembly considerations, 
etc.) were first imposed. This left three free parameters: 
ring electrode length, compensation electrode length, and 
end electrode length. Ring electrode length and end elec- 
trode length were chosen in order to both minimize Cq 
and achieve a large tunability with respect to C4, where 
tunability is described in Ref. [l7| : 



dC4 

Tunability = Vq— — . 

aVi 



(15) 



After determining the ring and end electrode lengths in 
this way, E2 was minimized with respect to the remain- 
ing parameter, the compensation electrode length, thereby 
"orthogonalizing" the geometry. Since E2 changes sign 
when scanning over electrode length, it was possible to 
choose a geometry that resulted in an arbitrarily small 
value for i?2- 

3. TAMUTRAP 

Performing the optimization as discussed above resulted 
in the geometry shown in Fig.|3l The trap radius is 90 mm, 
which is larger than in any existing Penning trap. The ring 
electrode length is 29.17 mm, the compensation electrode 
length is 71.36 mm, the end electrode length is 80.00 mm, 
and 0.50 mm gaps have been accounted for. This geome- 
try results in a length/radius ratio oi I/vq = 3.72. A good 
quadrupole field (C4 = -6.8 x lO^^, = 6.2 x 10"^) has 
been calculated to be achievable with compensation elec- 
trodes set to Vi ~ — O.373V0 (where Vq is the primary trap 
depth) . The analytic expansion of the electric field around 
the trap center up to Cs is shown in Table [21 Changing 
the voltage on the endcap electrodes (detectors) will adjust 
the predicted tuning (compensation) voltage; however, the 
trap will always remain "tunable" and "orthogonalized" 
since the contributions to the potential for each electrode 
are independent by superposition. 

4. Simulation 

The analytic solution for the proposed geometry was 
verified using SIMION 1^, an electric field and ion trajec- 
tory simulation program, in order to confirm the validity 
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Figure 3: The optimized trap geometry with SIMIDN generated elec- 
tric fields lines overlayed. 
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Table 2: Expansion coefBcients are compared for the optimized 
TAMUTRAP measurement trap when tuned (analytic and simulated) 
and two other existing Penning traps: TITAN (calculated analyt- 
ically as in Ref. [3), and LEBIT (simulated using SIMION as in 
Ref. [13). 



of the calculations. The results from SIMION are listed in 
Table [5] (along with the analytic solutions from |J3|) , and 
the resulting equipotential lines have been overlayed on 
the geometry cross section presented in Fig. [31 The val- 
ues output by SIMION agree with the analytic solutions for 
each Ck to a few parts in 10-3. The discrepancies between 
the analytic and simulated values can be accounted for by 
the inherent pixelation of the geometry as represented in 
SIMION and processing constraints, which have both been 
minimized as far as allowed by RAM and available com- 
puting time. Specifically, the optimized geometry has been 
represented in the simulation to the nearest 0.01 mm, and 
the voltages have been defined to 1 x 10-^ V. 

5. Comparison to Existing Traps 

The defining features of the TAMUTRAP geometry de- 
scribed here are the unique inner radius (90 mm), the 
small length/radius ratio (l/ro — 3.72), and the consider- 
ation made for detectors (short endcap electrodes). This 
new l/ro ratio is what allows TAMUTRAP to exhibit such an 
unprecedented radial size when compared to other cylin- 
drical Penning traps, while still measuring only 335 mm in 
overall length. The optimized ISOLTRAP geometry, for 
example, has an l/ro = 11-75 f21], and would therefore 
require over 1 m in length in order to maintain geometri- 
cal proportions (and, therefore, electric field shape) with 



a 90 mm inner radius. The compact geometry employed 
by TAMUTRAP allows for a structure with a very large radius 
to easily fit within the 1-m long bore of the 7T 210 ASR 
magnet. At the same time, the new analytic solution de- 
scribed in t|2. 31 retains the good quadrupole nature of the 
electric field displayed by other prominent Penning traps, 
which is required for other experiments such as precision 
mass measurements. 

Table [2] compares the analytic and simulated electric 
field expansion coefficients of TAMUTRAP to an analytic solu- 
tion reported by LEBIT and a simulated solution re- 
ported by TITAN 19] . The suppression of the anharmonic 
terms in the electric field generated by the geometry for 
the TAMUTRAP measurement Penning trap is comparable to 
that presented by these two prominent mass-measurement 
facilities, for which a very well-tuned harmonic electric 
field is critical [2^. With respect to the inherent field 
shape, the TAMUTRAP geometry should therefore be suitable 
for such precision mass measurements; however, it remains 
to be seen what affects the unprecedented electrode size 
and trapping volume necessitated by the primary program 
of performing /3 — v correlation measurements will have on 
the specific procedures required in these studies. In par- 
ticular, due to the enlarged geometry, the uniformity of 
the applied potentials across the electrodes is likely to be 
worse at TAMUTRAP than at existing facilities that employ 
smaller electrodes. It remains to be seen, however, whether 
or not these effects will ultimately limit mass resolution. 
In either case, they will not have a significant effect on 
the main program where the harmonic potential is not as 
stringent. 

6. Conclusion 

A novel, large-bore cylindrical Penning trap with a new 
length/radius ratio in the electrode structure has been de- 
scribed. The design work has been performed with the 
initial research program of measuring correlation parame- 
ters for T = 2 superallowed /3-delayed proton emitters in 
mind; however, careful attention has also been paid to cre- 
ating a facility with maximum suitability for a wide range 
of possible future nuclear physics experiments. In partic- 
ular, a line of sight to the trap center, an open geometry, 
and a well characterized quadrupole electric field are pre- 
dicted with the proposed design. The analytic solutions 
to the electric field were checked against simulation, and 
match to high precision. These values, which are com- 
parable to those presented by existing mass measurement 
facilities, suggest that TAMUTRAP should be capable of per- 
forming mass measurements in addition to its primary pro- 
gram. The aforementioned features, in conjunction with 
the unprecedented 90 mm free trap radius, will make the 
TAMUTRAP facility a unique tool for future nuclear physics 
research. 



5 



7. Acknowledgments 



The authors would hke to thank Guy Savard, Jason 
Clark, and Jens Billing for fruitful discussions and advice. 
This work was supported by the U.S. Department of En- 
ergy under Grant Numbers ER41747 and ER40773. 

References 

[1] N. Severijns, O. Naviliat-Cuncic, Annu. Rev. Nucl. Part. Sci. 61 
(2011) 23-46. 

[2] N. Severijns, M. Beck, O. Naviliat-Cuncic, Rev. Mod. Phys. 78 
(2006) 991-1040. 

[3] .1. A. Bchr, G. Gwinncr, J. Phys. G: Nucl. Part. Phys. 36 (3) 

(2009) 033101. 

[4] I. Towner, ,1. Hardy, Report on Progress in Physics 73 (4) (2010) 
046301. 

[5] N. Severijns, M. Tandecki, T. Phalct, I. Towner, Phys. Rev. C: 

Nucl. Phys. 78 (2002) 055501. 
[6] E. Adelberger, et al., Phys. Rev. Lett. 83 (7) (1999) 1299-1302. 
[7] J. Jackson, S. Treiman, H. Wyld, Nucl. Phys. 4 (1957) 206-212. 
[8] K. Blaum, Physics Reports 425 (1) (2006) 1-78. 
[9] M. Amoretti, et al.. Nature 419 (6906) (2002) 456-9. 
[10] V. Y. Kozlov, et al., Hyperfine Interact. 175 (2006) 15-22. 
[11] M. Mehlman, D. Melconian, P. D. Shidling, arXiv: 1208.4078 
(Aug. 2012). 

[12] Agilent Technologies, http://www.agilent.com (Aug. 2012). 

[13] C. Roux, et al., arXiv:1110.2920vl (Oct. 2011). 

[14] D. Beck, et al., Nucl. Instrum. Methods Phys. Res., Sect. B 

126 (96) (1997) 378-382. 
[15] G. Gabrielse, L. Haarsma, S. L. Rolston, J. Mass Spectrom. 88 

(1989) 319-332. 

[16] Wolfram Research, http://www.wolfram.com (Aug. 2012). 
[17] G. Gabrielse, Phys. Rev. A: At., Mol., Opt. Phys. 27 (5) (1983) 
2277. 

[18] Scientific Instrument Services, Inc., http://www.simion.com 
(Aug. 2012). 

[19] M. Brodeur, Ph.D. thesis. University of British Columbia 

(2010) . 

[20] R. Ringlc, et al., Nucl. Instrum. Methods Phys. Res., Sect. A 

3 (604) (2009) 536-547. 
[21] H. Raiinbault-Hartinann, et al., Nucl. Instrum. Methods Phys. 

Res., Sect. B 126 (1-4) (1997) 378-382. 
[22] G. Gabrielse, Int. J. Mass Spectrom. 279 (2-3) (2009) 107-112. 



6 



